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| Abstract 

> : 

Motivated by the observed forward-backward asymmetry in B — > K*£ + £~ decay, we 
perform a detailed analysis of this decay mode within a family non-universal Z' model. 
£Sj ■ With the related coupling Z' — sb constrained by B s — B s mixing, B — > ttK and B — > 

o ; 

\ X s [i + [i decays, we look for further constraint on the couplings Z' — from Aps{B — > 

A'* / u + /_f~) 0GcV 2^ 1J 2^2GeV 2 an d § e t numerically B^^ ~ O(10~ 2 ). Moreover, we find that 
the relations B^ < B^ and B^ + B^ < 0, with a small negative phase 4>si are 
crucial to moderate the discrepancy for Afb{B — > K*[i + [j>~) between the SM prediction 
and the experimental data. Numerically, comparing with the SM prediction, we find 
that Afb(B — > A*/x + ^ _ ) 0GcV 2^ (? 2 s g2GcV 2 com d be enhanced about 80% and 50% by Z' 
contribution at most in scenarios SI and S2, corresponding to the two fitted results of 4> s 
by UTfit collaboration, respectively. However, the results are still about 1.5a lower than 
the experimental measurement. 
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1 Introduction 



As is well-known, the electro-weak penguin decays b — > s£ + £~ appear only at the one-loop level 
in the Standard Model (SM), and are therefore very sensitive to possible new physics (NP) 
beyond the SM. Among many inclusive and exclusive processes based on the quark level b — > 
s£ + £~ transition, the exclusive B — > K*£ + £~ decays are of particular interest in this respect, 
as many observables in these decays, such as the branching ratio, the longitudinal polarization 
fraction, the forward-backward asymmetry Afb, and the isospin asymmetry, could be used to 
test the SM and to probe possible NP. 

In the literature, the exclusive B — > K*£ + £~ decays have been investigated in great detail 
by many authors, both in the SM and within various NP models [1, 2, 3, 4, 5, 6]. Among these 
observables of B — > K*£ + £~, the A FB is particularly useful. As discussed in Refs. [2, 3, 4], the 
zero of Afb is largely free from hadronic uncertainties in the SM and hence could be a powerful 
probe for various NP models. Recently this observable has been measured as a function of the 
dilepton invariant mass square q 2 = Mfi, by both BaBar [7, 8] and Belle [9, 10] collaborations. 
Their fitted Afb spectrum is generally higher than the SM expectation in all q 2 bins. Especially, 
the recent measurement from Belle collaboration [10] shows 

A FB {B K*£ + £-) 0GcVHqH2GcY , = 0.47^;1 ± 0.03, (1) 

favoring a positive value, whereas the sign of the SM prediction for Afb{B — > K*£ + £~) ~ —0.1 
at OGeV 2 ^ q 2 ^ 2GeV 2 is negative. Such a large discrepancy is hard to be moderated within 
the SM. Furthermore, the measurements prefer positive values for Afb{B — > K*£ + £~) in the 
whole q 2 region, indicating that there might be no zero crossing, which is apparently contrary 
to the SM prediction [2, 3, 4]. 

The measurements have motivated many recent investigations on the possible mismatch [11, 
12, 13]. In this paper, we revisit this decay mode within a family non-universal Z' model [14], 
which could be naturally derived in certain string constructions, E$ models and so on. In our 
previous paper [15], with the constraints from B s — B s mixing, B — > ttK and B — > X s/ u + /i~ 
decays, we have obtained an explicit picture for the Z 1 couplings B^ b and B^. Thus, it is 
of interest to see whether the discrepancy of Afb(B — > K*£ + £~) at ^ q 2 ^ 2GeV 2 between 
the SM prediction and the experimental data could be moderated with the constrained non- 
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universal Z' couplings. 

Our paper is organized as follows. In Sec. 2, we present a brief review of the SM theoretical 
framework for B — > K^*'£ + £~ decays. After brief introduction of the employed family non- 
universal Z 1 model in Sec. 3, we present our numerical analyses and discussions in Sec. 4. Our 
conclusions are summarized in Sec. 5. Appendices A and B include our theoretical inputs. 

2 The SM prediction 

Neglecting the doubly Cabibbo-suppressed contributions, the effective Hamiltonian governing 
b — > s£ + £~ transition is given by [5, 16] 

AC 10 

n cS = —j[V a VZ Yl CtMOiQi) , (2) 

where explicit expressions of Oi could be found in Ref. [5], and the Wilson coefficients C\ can 
be calculated perturbatively [3, 17, 18, 19], with the numerical results listed in Table 1. 



Table 1: The SM Wilson coefficients at the scale \i - 


= m b . 




C x {m b ) C 2 (m b ) C 3 (m b ) C 4 (m b ) C 5 (m b ) C 6 (m b ) Cf(m b ) Cf\ 


[m b ) - Y(q 2 ) 




-0.274 1.007 -0.004 0.076 0.000 0.001 -0.302 


4.094 


-4.193 


The effective coefficients Cf^ lQ in Table 1 are defined respectively 


as [20] 




^eff 4tt 1 4 20 80 
67 ~ o~ 7 ~ 3 3 ~ 9 4 ~ 3 5 ~ 9 ' 


3 ) 




Cg = — Cg + Y(q 2 ) , C^q = — C10 , 

a* a. 




(3) 



where Y(q 2 ) denotes the matrix element of four-quark operators and is known from the liter- 
ature [4, 5, 21, 22]. We have neglected the long-distance contribution mainly due to J /\& and 
ty' in the decay chain B — > lf(*)\]/W — y K { ^ ) t r t~ , which could be vetoed experimentally [8, 10]. 
For recent detailed discussion of such resonance effects, we refer to Ref. [23]. 

Although there are quite a lot of interesting observables in B — > K^*H + £~ decay, in this 
paper we shall focus only on the dilepton invariant mass spectrum and the forward-backward 
asymmetry. 
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Adopting the same convention and notation as [4], the dilepton invariant mass spectrum 
for B -> K(*H + £- decay is given respectively as [4, 24] 



dr 



K 



ds 



G% a 2 m 5 



us 



C'\ 2 4mj(2 + 2m 2 K - s) 



+Re(C'D'*)8mj(l - m 2 K ) + \D'\ 2 Amjs 



(4) 



dr 



K " 



G% a 2 ml 



B 



ds 



2 10 7T 5 
1 

4m 2 K , 




A\ 2 AW m 2 p. . w(s) 



Am 2 K , 



2m 2 K * 



\B\ 2 (\ 
\C\ 2 {\ 



us 



8m^(s + 2m^)) 



| F |2 (A _!H£L + 8m^(s 



|Gf A- 



us 



4m 2 (2 + 2m 



Re(BC*)(A 



us 



-)(l-m|*- S ) 



+ Re(FG*)((A- 



u s 



■)(1 -fflx. - s) +4m^A) 



-2- 



■A [Re(FH*) - Re(GH*)(l - m 2 K ,)] 



m c 



■s\\H\ 2 



m 



(5) 



K" 



Here the auxiliary functions with the explicit expressions given in [4], are combinations 

of the effective Wilson coefficients Eq. (3) and the B — > K^*' transition form factors, which are 
calculated with light-cone QCD sum rule approach in Ref. [25]. 

The differential forward-backward asymmetry for B — > K*£ + £~ decay is defined as [4] 



dw4.pg 



G\a 2 m\ 



llb B I T T T |2 ~ ~/-\2 



X 



ReiC^C^VA, + !^Re(C , 7 Gff C 1 c ff *) (VT 2 (1 - m K .) + A 1 T 1 (1 + m K « 



.(6) 



It is noted that, although C^q is real in the SM, it could become complex after including the 
Z' contributions. Note that the Afb for B — > K£ + £~ decay vanish both in the SM and within 
the Z' model considered in the present paper, since neither further operators nor higher-order 
corrections are included in our discussion [4, 26]. From the experimental point of view, the 
normalized forward-backward asymmetry is more useful, which is defined as [4] 

d^FB 



ds 



d*<4.FB , dr 
ds ds 



(7) 



3 Family non-universal Z' couplings and their effects 

A family non-universal Z' model, which has been formulated in detail by Langacker and 
Pliimacher [14], can lead to FCNC processes even at tree level due to the non-diagonal chiral 
coupling matrix. Under the assumption that the couplings of right-handed quark flavors with 
Z 1 boson are diagonal, the Z' part of the effective Hamiltonian for b — > sl + l~ can be written 
as [15, 27] 



2C 

H z e 'Ab^sl + r) = £v tb V? s -^^{sb) v _Mv-A-^f^{sb) v _ A {ll)v +A +h.c, (8) 



nL pi tdL tdR 



V2 



which could also be reformulated as 



V tb V t * s * Vtb v ts 



AC 

H? ff (b sl+r) = ^V tb V* [AC 9 9 + AC[ O w ] + h.c. , (9) 



V2 



with 



e v ts v tb 



sb I r>L i tjR\ 



*C 10 (M W ) = 9 4J^r( B u- B $)> ( 10 ) 

e v ts V tb 

where B^ b and B^' R denote the effective chiral Z' couplings to quarks and leptons. 

With the above expressions, the Z' contributions can be represented as modifications of the 
Wilson coefficient of the corresponding semileptonic operators, i.e., C 9 10 (Mw) = Ci%{Mw) + 
AC' gl0 (Mw)- The running from Mw scale down to m b is the same as the SM ones[5, 16]. 
Numerically, with the central value of the inputs, we get 

C' 9 {m b ) = 0.0682 - 28.82 J3L(b£ + B*), (11) 

v ts v tb 

C[ (m b ) = -0.0695 + 28.82-^(^-5*), (12) 

vts V tb 

To include Z' effects, one just needs to make the replacements 

Att 

Cf _> Of = ^C' 9 + Y(q 2 ) , 
Att 

Cf -+ Cf =-C' w . (13) 

in the formalisms relevant to B — > K^l^l" listed in section 2. With the formulae collected 
above, we shall proceed to present our numerical analyses and discussions in the next section. 



4 Numerical analyses and discussions 

The considered Z' contributions to B — > Kfi + fi~ and B — y K*fi + fi~ decays involve three 
couplings: B^ b , B^ and Bj£ . In our previous papers [15], we have combined B s — B s mixing, 
B — y 71 and pK, B as well as B s — y decays to constrain these Z' 

couplings and their possible phase. Our combined results are re-tabulated in Table 2 and re- 
displayed in Fig. 1 (pink region). The two solutions SI and S2 for \B^ b \ and <f>g correspond to the 
two fitted results for the new physics parameter (f>B s performed by the UTfit collaboration [28]. 
It is natural to question whether the constrained parameter space could account for the Ap B 
measured recently by the Belle collaboration [10]. 

The Afb(s) spectrum for B —y K* decay measured by Belle collaboration tends to 
be shifted toward the positive side in all six q 2 bins, indicating that there might be no zero 
crossing. However, a zero crossing in Afb{s), whose position is well-determined and free from 
hadronic uncertainties at the leading order in a s , is well predicted in the SM [2, 3, 4]. Especially 
in the bin OGeV 2 < q 2 ^ 2GeV 2 , the sign of the SM prediction is negative, being different from 
the experimental measurement, which favors a positive value on the other hand. 

Comparing the two terms in the square bracket in Eq. (6), one can see that, at low q 2 region 
the first term Re(C$ eS C°Q *) is suppressed by one power of q 2 /m 2 relative to the second one, and 
A FB {B — y K*p + fi~) is therefore dominated by the second term Re(C 7 eG CiQ *). Thus, at low q 2 
region A FB (B — y K* fi + p:~) could be significantly changed by flipping the sign of Re(Cy e *) 
from a positive in the SM (due to C^ s < and Cfg < 0) to a negative one. On the other hand, 
in order to keep Afb(B — y K*p + fi~) positive in the high q 2 region, the sign of Re(Cg e CJo*) 
should be maintained negative as predicted in the SM. 

With 0f 72° (SI) and <p L s 82° (S2) obtained in Ref. [15], and keeping in mind that 

the CKM element Vt s is negative, one can easily find from Eq. (12) that the sign of Re(C7 e C°q*) 
could be flipped if B(( < B^. In order to see the Z' effect on Re(Cg S C^Q *) explicitly, with Y(q 2 ) 
excluded, we can rewrite it as 

Re(CfC$*) = Re(Cf ) Re(C 1 c ff *) + Jm(Of) Jm(0$) , 

~ Re(Cf ) Re(Cf *) + ^ * Im(A^) Im( AC( ) , (14) 

where the fact that both Cg G and C^q are real in the SM has been used in the second line. 
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Table 2: Columns 2-5 present the values for B^ b with the constraints from B s — B s mixing and 
B — > ttK*, pK decays and Bz^ 1 with the constraints from B — > X s fi + fi~ decay [15]. Columns 
6-7 are the re-fitted results after including the constraints from B — > K*fi + p~ decay. 





|£4l(xl0- 3 ) 




£^(xl0" 2 ) 


<(xl0" 2 ) 


^(xlO- 2 ) 


^(xlO- 2 ) 


SI 


1.09 ±0.22 


-72 ±7 


-2.7 ±2.5 


0.61 ±2.4 


-4.75 ± 2.44 


1.97 ±2.24 


S2 


2.20 ±0.15 


-82 ±4 


-0.59 ±0.93 


0.19 ±0.88 


-1.83 ±0.82 


0.68 ±0.85 



To keep the sign of Re((7g ff )Re(C^o*) negative, one can derive the relation Bfc + Bjj < from 
Eq. (11). At the same time, with the obtained relations, B^ + B^ < and B[j < Bjf, the term 
Im(AC , g)Im(AC , { ) is automatically negative, and hence the sign of Re(C 9 e// C 1 e f*) is indeed 
maintained to be negative. It is interesting to note that the allowed parameter space in Table 2 
constrained by B — > X s fi + fi~ and B s — > p + fJ>~ decays [15] satisfy the relations, Bfj < Bff and 
Bf( + Bft < 0. It is however unclear whether the parameter space could bridge the discrepancy 
of Afb{B — > K*fi + fi~) between the SM prediction and the experimental data. In the following 
numerical evaluation, we shall perform a fit combining the constraints from B — > X s fi + fi~ and 
A FB (B — > -ft"*/i + /i~)oGcV 2 s£g 2 s£2GeV 2 ! leaving the other observables for B — > p + /i~ decays as 
our predictions within such a Z' model. 

For consistence, we take the same simplifications for the family non-universal Z' couplings as 
Ref. [15]. Our fit is performed with the experimental data on Afb(B — > K* fi + fi~) 0Gc y2^ q2 ^ 2 GeV 2 
varying randomly within its 1.7cr (~ 90% C.L.) error bar, while the theoretical uncertainties are 
obtained by varying the input parameters within their respective regions specified in Appendices 

A and B. Our re-fitted numerical results for B,,:, are listed in the columns 6-7 in Table 2, and 

mm ' 

the corresponding allowed regions are shown in Fig. 1 (green region). As illustrated in Fig. 1, 
after including the constraint from A FB (B — > K* /i + /i~)oGcV 2 ^ 2 s£2GeV 2 5 the survived parameter 
spaces of Z' couplings (green region) is further reduced. Comparing with the constraint from 
B —)• X s fi + fi~ decays (pink region) only, one can see that the regions with B^ > B^ are 
excluded, which confirms our naive analysis that the relation Bjf < B^ is needed to alleviate 
the observed discrepancy for A FB (B — > K*fi + fi~) QGc y2^ q 2^ 2 Gcy 2 - Numerically we also find that 
B^ < and \B^ \ > \B^ \, which means that the other condition B^ + B^ < is also kept. 
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Figure 1: The allowed regions for the parameters B^lf 1 . The pink regions are allowed by the 
constraints from B — > X s/ u + /i~ decay [15]. The green ones correspond to the final parameter 
space with the constraint from Aps(B — > K* /i + /W~) 0GcV 2^ g 2^2GeV 2 a ^ so included. The shaded 
region is derived from our qualitative analysis with the conditions Bh + B^ < and B^ < B^. 



With the constrained Z' couplings and taking q 2 = lGeV 2 , we show in Fig. 2 the dependence 
of A FB (B — > K*fi + fi~) on \Bg b \, 0f and B^jf 1 . From Fig. 2 (a) and (b), one can see that 
Afb{B — > K*fi + fi~) at low q 2 bin could be enhanced to the experimental side by the Z' 
contribution with a large negative Bjj and/or a positive Bff . At the same time, from Fig. 2 (c) 
and (d) one can see that, with the constrained \B^ b \, a smaller phase |0f | is more helpful to 
enhance the A F b(B — > K*n + n~) at low q 2 region. As a result, due to the fact that [(/^[(Sl) < 
\(pg\(S2), the solution SI is preferable to S2. Moreover, a larger \B^ b \, which enlarges the Z' 
contribution, is also helpful for reconciling the A FB {B — > K*fi + fi~) problem at low q 2 . 

With the preferred choice \B^\ = 1.1 (2.2) x 10~ 3 , <f% = -65° (-78°), B L m = -7.2 (-2.7) x 
10 -2 , B^ = 4.2 (1.5) x 10~ 2 , and the central values of the other inputs, we get 

A FB (B -> A> + /i-) 0GeV2 ^ 2GeV2 = -0.02 (-0.05) (15) 

in scenario SI (S2). Compared with the SM prediction ~ —0.10, this observable could be 
enhanced about 80% (50%) by the Z' contribution, implying that the scenario SI with <pB s = 
—20.3° ± 5.3° fitted by UTfit collaboration [28] is favored by these decays. However, since 
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Figure 2: The dependence of Afb(B 
comparison, the SM prediction is shown as dashed 



on \BL\, 
lines. 



j> L s and BfrK with q 2 = lGeV 2 . For 



B — > X s fi + fi decay has already put strong constraint on the strength of the Z' couplings, the 
result for A FB (B -> K* 

A^A* )oGeV 2 <g 2 ^2GcV 2 ^ s still negative in SI (S2), being 1.5 (1.6)<r lower 
than the data QA7t° 2 £, [10]. Such a situation could also be seen from Fig. 3, where the effects 
of Z' contribution induced by B^ R on dB{B ->■ K^n + yr)/ds and A FB (B ->■ K*fi + n~) are 
displayed. 

With the bounded Z' couplings listed in Table 2, our predictions for B(B — > K^*> /i + 
and Afb{B — > K*fi + fi~), both within the SM and in the Z' model, are given in Tables 3 
and 4. We find that most of the observables agree with the experimental data within errors. 
However, beside the discrepancy of Afb{B — > K*fi + fi~) in the low q 2 bin, a similar problem 
is also observed in the high q 2 bin (q 2 ^ 16GeV 2 ). Within the allowed parameter space shown 
in Fig. 1, the discrepancy for Afb(B —> K* ji + yU~) 9 2^ 16GeV 2 between the SM prediction and 
the experimental data is still difficult to be reconciled. In fact, if the Z' correction could not 
give a significant contribution after totally counteracting the SM contributions, the Afb{B — > 
K*fi + fi~) problems in both the low and the high q 2 bins will persist within such a family 
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Table 3: Predictions for B(B — y Kfi + fi ) within the SM and the Z' model. 



q 2 (GeV 2 ) 


Exp. [10] 


SM 


SI 


S2 


[0,25] 


4.5 ± 0.4 [30] 


5.57 ±0.51 


5.74 ±0.84 


6.16 ±0.63 


[1,6] 


1-36^ 


1.62 ±0.21 


1.64 ±0.29 


1.77 ±0.24 


[0,2] 


0.81jft? 


0.62 ±0.09 


0.61 ±0.12 


0.66 ±0.10 


^ 16 


o-98±8:S 


0.84 ±0.03 


0.86 ±0.11 


0.92 ±0.07 



non-universal Z' model. 

As a final comment, abandoning the constraints from B — y X s fi + fi~ decay, we pursue the 
required strength of Z' couplings in order to make the sign of Ke(C^ s C^*) flipped, while leaving 
Re (CfC$*) unchanged. The conditions Re(CfC$*) < and Re(CfC$f*) < are equivalent 
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Table 4: Predictions for B{B ->■ K*fi + fi~) and A FB (S -> K*n + fi~) within the SM and the Z' 
model. 





q 2 (GeV 2 ) 


Exp. [10] 


SM 


SI 


S2 




[0, 25] 


10.8±H [3°] 


11.3 ± 0.2 


11.8 ± 1.4 


12.6 ± 0.9 




[1,6] 


1 4Q+0.47 


2.35 ±0.05 


2.45 ±0.31 


2.61 ±0.19 




[0,2] 




1.23 ±0.06 


1.24 ±0.12 


1.26 ±0.08 




^ 16 


2 04 +0 ' 32 


1.39 ±0.03 


1.45 ±0.18 


1.55 ±0.12 




[0,25] 




0.27 ±0.01 


0.25 ±0.05 


0.26 ±0.03 




[1,6] 


26+ 28 


0.07 ±0.01 


0.09 ±0.03 


0.09 ±0.02 




[0,2] 


n 47+ - 26 

u -^' -0.32 


-0.10 ±0.01 


-0.05 ±0.02 


-0.07 ±0.01 




^ 16 


66+ 012 


0.34 ±0.01 


0.30 ±0.05 


0.32 ±0.03 



to the following demands: 

Re(C%) > 0, Re(Cg) < Im(AC 9 ) > . (16) 

From Eqs. (11) and (12), and with the values of \B^ b \ and listed in Table 2, we get 

< -0.271 , 0.273 ± B^ < B% < -0.268 - B^ , (SI) 
B L m < -0.287, 0.280 ± B L m < B% < -0.275 - B^ , (S2) (17) 

Taking B^ = —0.3 and B^ = 0.01 in both SI and S2, which satisfy Eq. (17) and are similar to 
the results given by Ref. [11], as shown in Fig. 3 (c) Afb(B — > K*fi + fi~) could be significantly 
enhanced to the experimental level. Unfortunately, such a region given by Eq. (17) is excluded 
by the constraint from B — > X s ji + ji~ decay. 

5 Conclusion 

In conclusion, motivated by the large discrepancy for Afb{B — > K*fi + fi~) in the low q 2 region, 
we have studied a family non-universal Z' model to pursue possible solution. With the con- 
strained coupling Z' — sb from B s — B s mixing and B — > ttK decays [15], we focus on the allowed 
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regions for the couplings Z' — fi + fi Bjjjf', which have already been strongly constrained by 
BiB — > X s fi + fi~) decay [15] at high, low and full q 2 regions, with the further constraint from 
A FB (B -)■ K* H + /i )oGeV 2 ^ij 2 ^2GcV 2 • Within the allowed Z' couplings, we have investigated the 
effect of such a Z' model on the observables of B — > /i + pT decays. Our main conclusions 
are summarized as: 

• To account for the experimental data on A F b{B — > K*fi + fi~), naively we get two 
interesting relations, B^ < B^ and B^ + Bj£ < 0. Furthermore, a larger \B^ b \ 
and a smaller |0f;| (negative) are crucial to moderate the discrepancy of A FB {B — > 
K* /i + /i~)oGeV 2 s£(} 2 s£2GeV 2 between the SM prediction and experimental data. Thus, sce- 
nario SI is preferable to S2. 

• A FB {B -> K* 

A* "'"A* )oGeV 2 s£<j 2 s£2GeV 2 puts a strong constraint on the Z' couplings, B^^ . In- 
cluding the constraints from B — >■ X s fi + fi~ decay and A FB (B — > K*{i + fi~) 0GcV 2^ q 2^ 2Ge y2, 
we get B^ + ^ ~ -5(-2) x 10" 2 and B* +fi _ ~ 2(1) x 10" 2 in SI (S2). 

• Due to the severe constraints from B — > X s fi + fi~ on the strength of the Z' contri- 
bution, the A F b{B — y K*fi + fi~) problems in both low and high q 2 regions still per- 
sist. After including the Z' contribution, compared to the SM prediction, A F b{B — > 
-Ar*/i + /i~) 0GeV 2 s g (? 2 s g2GeV 2 could be enhanced by an amount of about 80% (50%), which is 
still about 1.5a (1.6<r) lower than the experimental data in SI (S2). 

Within such a family non-universal Z' model, although involving the same Z' — sb coupling, 
these different processes also depend on different diagonal Z' couplings. For example, B — > 
X s /j + [x~ and B — > K^n + /j~ depend on lepton diagonal coupling Z' — while B — > Kir 

decays on quark diagonal couplings Z' — uu and Z' — dd. So, if one observable from the list of 
current anomalies (A F b at low dilepton mass, the B s mixing phase, the Ktt puzzle) becomes SM- 
like, it will only affect the respective coupling, while leaving the others unchanged. For example, 
if the B s mixing phase becomes SM-like, then through adjusting the diagonal lepton and quark 
couplings, we can still find suitable parameter spaces to account for the other anomalies. If 
the Ktt puzzle and/or the A F b at low dilepton mass become SM-like, it will only give more 
severe constraints on the diagonal lepton and quark couplings. To further constrain the model 
parameter spaces, it is therefore necessary to combine all of these processes at the same time 
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and perform a global analysis, which is however beyond the scope of this paper but will be 
addressed in a forthcoming publication. With the upcoming LHC-b and proposed super-B 
experiments, the data on these processes is expected to be more precise [31], which will then 
severely shrink or totally excluded the model. 

Note added: During our work on the way, we note that a recent paper [11] also pursues possible 
solutions within a family non-universal Z' model. In order to enhance the Z' contribution 
to the real part of Cf^Cf^*, they have assumed that 0f = 0, which is obviously unsuitable 
due to the fact that a nonzero 0f is needed to resolve the "ttK puzzle" and the U B S — B s 
problem" [15, 27, 29]. Furthermore, their result Bf[ ~ Bjj ~ (^(lO" 1 ) is also excluded by the 
constraint from B — > X s p + p~ decays, in which Bfc ~ Bff ~ (9(10~ 2 ) [15]. 
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Appendix A: Theoretical input parameters 

For the CKM matrix elements, we adopt the fitting results from the UTfit collaboration [28, 32] 

p = 0.154 ± 0.022 (0.177 ± 0.044), 
fj = 0.342 ± 0.014 (0.360 ± 0.031), 
\V U /V ts \ = 0.209 ± 0.0075 (0.206 ± 0.012), 

\V cb \ = (4.13 ± 0.05) x 10~ 2 ((4.12 ± 0.05) x 10~ 2 ), (18) 

with p — p (1 — 4r) and fj — r\ (1 — 4r). The values given in the brackets are the CKM parameters 
in presence of generic NP, and are used in our calculation when the Z' contributions are included. 
As for the quark masses, we take [33, 34] 

m u = nid = m s = 0, m c = 1.6ljT 01 | GeV, 

m b = 4.79±g;5i GeV, m t = 172.4 ± 1.22 GeV. (19) 
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Appendix B: Transition form factors from light-cone QCD 
sum rule 



In order to calculate the B — > K^*'£ + £~ decay amplitude, we have to evaluate the B — > 
matrix elements of quark bilinear currents. They can be expressed in terms of ten form factors, 
which depend on the momentum transfer q 2 between the B and the mesons (q = p—k) [25]: 



(K(k)\s lt ,(l- l5 )b\B{p)) 

(K(k)\sa, u q»(l + l5 )b\B(p)) 
with / + (0) = / (0), 

(K*(k)\s lfl (l- l5 )b\B(p)) 



U{q 2 ) (2p-ql 



in 



B 



m 



K 



2 _2 



m. 



(2p-q)^q -q^(m B -m K 



+ 



f T (q 2 



m 



K 



(21) 



-ie* (m B + rriK^A^q 2 ) + i(2p - q)^(e* ■ q) 



A 2 (q 2 



mB + m K * 



+iq^(e* ■ q) 



A 3 (q 2 ) - A (q 2 ) 



2V{q 2 



ituupat P K> 

tub + rriK* 

with A 3 (q 2 ) = H^^ Al (q 2 ) - «i 2 ( 9 2 ) and 4,(0) = A 3 (0), 



(22) 



(K*(k)\sa^(l + l5 )b\B(p)} 



te^e^pPk^T^q 2 



+T 2 (q 2 
+Uq 2 W-q) 



e^m B - m K * 



2 2 

m B — m K * 



(e*-g) {2p-q) fl 



{2p~q)» 



(23) 



with Ti(0) = T 2 (0). e M is the polarization vector of the K* meson. The physical range in s = q 2 
extends from s min = to s max = {m B — m K (*)) 2 . 

These transition form factors have been updated recently within the light-cone QCD sum 
rule approach [25]. For the q 2 dependence of the form factors, they can be parameterized in 
terms of simple formulae with two or three parameters. The form factors V, Aq and T\ are 
parameterized by 

ri r 2 



F{s) 



+ 



(24) 



l-s/m 2 R l-s/m 2 t 

For the form factors A 2 , T 3 , / + and fr, it is more appropriate to expand to the second order 
around the pole, yielding 

n r 2 



F(s) 



1 



s/m 2 (1 
14 



s/m)' 



(25) 



Table 5: Fit parameters for B — > transition form factors [25]. 





F(0) 


' 1 


2 


Vet 

1 2 


9 

Tfl 7, 

'"fit 






fB^K 
J + 


0.331 


0.162 


5.41 2 


0.173 




Ea 


(25) 


Jo 


0.331 






0.330 


37.46 


Ea 


(26) 


JT 


0.358 


0.161 


5.41 2 


0.198 




Ea 


(25) 


yB^K* 

X 


0.411 


923 


5 32 2 


—0 51 1 


49 40 


Ea 


(24) 


A B^K* 
u 


0.374 


1.364 


5.28 2 


-0.990 


36.78 


Eq. 


(24) 


A B->K* 


0.292 






0.290 


40.38 


Eq. 


(26) 


aB^K* 
A 2 


0.259 


-0.084 




0.342 


52.00 


Eq. 


(25) 


rpB^K* 


0.333 


0.823 


5.32 2 


-0.491 


46.31 


Eq. 


(24) 


rpB^K* 

1 2 


0.333 






0.333 


41.41 


Eq. 


(26) 


rfB^K* 

J 3 


0.333 


-0.036 




0.368 


48.10 


Eq. 


(25) 



where m = m%t for A-i and T3, and m = m,R for / + and fx- The fit formula for Ai, T2 and /o is 

FW = i^k- (26) 

The form factor T3 can be obtained through the relation Ts(s) = 1 ~™ K * [T^s) — ^(s)] . All the 
relevant fitting parameters for these form factors are taken from Ref. [25] and are recollected 
in Table 5. 
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